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Abstract. The recent global fit of short baseline neutrino oscillation data favors the presence 
of one (or more) sterile neutrino state which leads to new mass splitting Am 2 ~ 1 eV 2 . We 
' consider the effect of this new states on the evolution of neutrinos from the dark matter anni- 

hilation inside the Sun. We show that neutrinos with energy E v > 100 GeV undergo resonant 
active-sterile oscillation which depletes the flux of neutrinos arriving at the Earth. As an 
CN \ example of this effect, we present the oscillation probabilities for the case of monochromatic 

neutrinos from the direct annihilation of dark matter particles to neutrinos and the depletion 
due to the presence of sterile neutrinos. 
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1 Introduction 

Well-established evidences ranging from galaxy to clusters of galaxies and cosmological scales 
are in favor of the existence of a huge amount of matter in the universe, the so-called Dark 
Matter (DM) [1]. In spite of the compelling evidences on the existence of DM, almost 
nothing is known about the nature of DM and especially the fundamental particle(s) that 
can comprise it. Basically two approaches are used in experiments aiming at the detection of 
DM: direct and indirect methods. The direct detection method is based on measurement of 
the scattering of particles comprising DM (hereafter x) off the nuclei of ordinary matter. The 
current upper bounds, from experiments aiming at direct detection, on the spin-independent 
a?} N and spin-dependent cr^° cross sections are respectively ~ 10~ 8 pb [2] and ~ 10 -2 pb 
[3, 4] for DM mass m x ~ 100 GeV. 

The indirect method is based on the detection of annihilation products of DM particles 
including gamma rays, charged particles and neutrinos. To attain a substantial rate of XX 
(where x = X m the case of a Majorana DM particle) annihilation, the density of DM should 
be sizable which is the case in the center of the galaxies and astrophysical objects, and the 
closest one to us is the Sun which we consider in this paper. The DM particles which scatter 
elastically off the nuclei in Sun lose a part of their momentum and if their velocity drop 
below the escape velocity in the Sun, they trap gravitationally and accumulate in the center 
of the Sun (see [5] and references therein). Accumulation of the DM particles trapped by 
the gravitational potential well of the Sun increases the number density n x and since the 
annihilation rate T oc n 2 , large values of n x can result in sizable annihilation rate. 

Annihilation of DM particles inside the Sun can produce neutrinos in two ways: 1) As 
the secondary product of the annihilation modes XX ~^ QQ, fermion pairs, gauge bosons or 
gauge bosons and Higgs particle. In this case the neutrino beam has a continuous energy 
spectrum that its shape depends on the annihilation mode; 2) As the direct annihilation of 
the DM particles XX ~ > v v- I n this case the neutrino beam is monochromatic with energy 
equal to the DM particle's mass. It should be noted that although the branching ratio of the 
annihilation mode XX ~^ u ^ is predicted to be small in some models (especially in MSSM 
where the annihilation cross section of the lightest neutralino to neutrinos is suppressed by 
the factor m 2 /m 2 ), it is possible to build models that this suppression is circumvented (see 
[6, 7] for such models). It can be shown that if the DM-nuclei cross section saturates the 
upper bound from the direct method experiments, the number density n x is large enough to 
produce detectable flux of neutrinos in neutrino telescopes such as IceCube [8] . 

The recent global fit of short baseline neutrino oscillation data, including LSND [9], 
MiniBooNE [10] and reactor anomaly [11], shows indication in favor of active-sterile oscil- 
lation with mass splitting Am 2 ~ leV 2 [12]. In this paper we consider the effect of this 
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new mass splitting on the evolution of neutrino beam from DM annihilation inside the Sun. 
We show that the presence of a mass splitting Am 2 ~ 1 eV 2 would lead to resonant flavor 
conversion for neutrinos with energy E u > 100 GeV; similar to the resonant flavor conversion 
of thermonuclear solar neutrinos with E u ~ MeV and Am^ ~ 8 x 10 _5 eV 2 . We show that 
due to this resonant flavor conversion, the electron neutrinos and muon/tau anti-neutrinos 
produced in the annihilation of DM particles with mass m x > 100 GeV are mostly sterile 
when they leave the Sun; which results in a reduction of the expected number of events in 
neutrino telescopes such as IceCube [13]. 

As an example of the consequences of this resonant flavor conversion, we consider the 
case of monochromatic neutrinos from the annihilation XX ~ * y v. It was proposed in [8, 14] 
that in the case of monochromatic neutrinos the number of events at neutrino telescopes on 
the Earth (such as IceCube) would show a seasonal variation due to the change in the Earth- 
Sun distance during the year. We calculate the flavor oscillation probabilities from Sun to 
the IceCube site for monochromatic neutrinos and the expected deviations in the seasonal 
variation due to the presence of Am 2 ~ leV 2 . The analysis for the case of secondary 
neutrinos with continuous energy spectrum from the annihilations XX ~~ QQ> W + W~ , . . . will 
appear in [15]. 

This paper is organized as follows. In section 2 we discuss the flavor oscillation of 
neutrinos inside the Sun and the conditions that a resonant flavor oscillation would take 
place. In section 3 we calculate the oscillation probabilities at the Earth and we discuss 
the deviations in the seasonal variation resulting from the active-sterile oscillation. The 
conclusions are given in section 4. 



2 Oscillation inside the Sun 

The presence of new sterile neutrino states with mass splitting Am 2 ~ 1 eV 2 favored by the 
global analysis of short baseline neutrino oscillation data would change the evolution pattern 
of neutrinos inside the Sun. Especially, the matter effect of the Sun's medium would induce a 
MSW resonance [16] which can drastically change the flavor content of neutrino beam at the 
surface of the Sun. In this paper we assume one new mass eigenstate v± (the so-called 3+1 
scenario [17]) with mass splittings Am 2 x ~ Am| 2 — ^ m %, ~ leV 2 , where Am^ = m 2 — m 2 . 
To calculate the flavor oscillation probabilities we use the following convention for the mixing 
matrix (assuming all the CP-violating phases equal to zero): 

U = R34R24R14R23R13R12 (2.1) 

where Rij (with i,j = 1,2,3,4 and i < j) denotes the 4-dimensional rotation matrix in i-j 
plane with the rotation angle 6ij (we use the notation = sin 9ij and Cij = cos 6ij ) . For the 
conventional mixing angles #12, #23 and #13 we use the values from [18]. The current best-fit 
values for #14 and #24 from the global analysis of short baseline neutrino oscillation data [12] 
are: s 2 4 = 0.03 and s 2 ,^ = 0.012. For the mixing angle #34 the strongest limits comes from 
the MINOS experiment [19]: s| 4 < 0.09. 

The evolution equations of the flavor neutrino states inside the Sun are 



dr 



'UM 2 U^\ + V(r) 



2E V 



v p , (2.2) 



where r is the distance from Sun's center; a, (3 = e, fi, r, s; and 

M 2 = diag (0, Am 2 ,!, Am^, Am^) 
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where Am^ ~ 8 x 10~ 5 eV 2 and Am^ ~2x 10~ 3 eV 2 [18]; and 

V(r) = ^G F diag(Af e (r),0,0,^ n (r)/2). 

AT e and N n are respectively the electron and neutron number density profiles of the Sun. 

In the two flavors approximation, the MSW [16] resonance energy for the active-sterile 
oscillation is given by 

, „ /Am?,\ /cos 29 \ f N A cm~ 3 \ , . 

E -= 6 - 55x10 GeV (-n#)(— )(^r)- (2 - 3) 



where Na is the Avogadro's number and 9 is #14, $24 and #34 for e — s, \i — s and r — s 
resonance, respectively. However, due to the hierarchical mass scheme in 3+1 scenario and 
the experimental bounds on the active-sterile mixing angles, the two flavors approximation 
for the resonance energy is quite accurate. It should be emphasized that although we use 
eq. (2.3) for approximating the resonance energy, we calculate the probabilities in the realistic 
four flavors scheme numerically. For e — s resonance: N e ff = N e — N n /2 and for (/i/V) — s 
resonance: N e ff = N n /2. Also, we assume normal hierarchy for the active-sterile mass 
splitting; i.e., Am^ > 0. It should be noticed that the e — s resonance occurs for neutrinos 
and (///t) — s resonance occurs for the anti-neutrinos. The minimum value of the resonance 
energy, for fixed values of Am^ and cos 29, corresponds to the maximum value of N e jf which 
is at the center of Sun. According to Standard Solar Model [20], at the center of Sun we have 



Max[(AT e// r s ] =(N e ^ 



2 



N A 

76.9 



cm^ 

center 



Max 



(Neff) 



N A 

26.9 



center Cm3 



Thus, it is easy to see from eq. (2.3) that the evidences from short baseline neutrino exper- 
iments on the existence of a new sterile state with Am^ ~ 1 eV 2 and small mixing angles 
cos2#j4 ~ 1 {i = 1,2,3), would lead to resonant flavor conversion v e — > v s for electron neu- 
trinos with energy E v > 85 GeV and P^/ T — > v s for muon and tau anti-neutrinos with energy 
£„> 240GeV. 

To illustrate the resonance behavior we plot in figure 1 the effective mass-squared differ- 
ence in matter (Am^M with respect to the radial distance from the center of Sun. In this 
plot we assumed: s 2 4 = 0.03 and Am^ = leV 2 . The red (dashed) and blue (dotted) curves 
correspond respectively to neutrinos with energy E v = 200 GeV and 600 GeV. The minima 
in these curves show the resonance point in the trajectory of electron neutrinos toward the 
surface of Sun. At the same plot of figure 1, the black (solid) curve shows the number density 
of N e — N n /2 (with the scale at the right vertical axis) which plays role in the e — s resonant 
flavor oscillation. As can be seen the resonance occur at the distance from the Sun's center 
where the condition in eq. (2.3) is satisfied. The DM particles trapped inside the Sun, would 
virialize with the Sun's medium. Comparing the gravitational potential energy of the Sun's 
matter with the average thermal energy of DM it can be shown that DM particles reside in 
the center of Sun in a radius tdm % Q-OIRq. On the other hand, according to the Standard 
Solar Model, the matter density in the Sun is almost constant within the radius 0.05Rq, 
which means that all the neutrinos produced from XX annihilation with energy higher than 
the resonance energy in eq. (2.3) would undergo MSW active-sterile resonance. 
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Figure 1. The effective mass-squared difference in matter (Am 41 )M (left vertical axis) and N e — N n /2 
(right vertical axis) with respect to the distance from Sun's center. The red (dashed) and blue (dotted) 
curves correspond respectively to neutrinos with energy E v = 200 GcV and 600 GeV at the center of 
Sun and the black (solid) curve shows the number density. In this plot we assumed: sin 2 814 = 0.03 
and Ami, = lcV 2 . 




E v (GeV) 



Figure 2. The probabilities P{vf ^ 0surf ) (red curve) and P{vf vf smt ) (blue curve) with 
respect to the neutrino energy. In this plot we assumed (sin 2 #14 = 0.03, 9 2 a = 0, 6*34 = 0) and 
(6»i4 = 0,sin 2 6»24 = 0.01,6*34 = 0) for -> ^ 0surf and vf zy 0surf , respectively; and for the both 
cases Am^ = 1 eV 2 . 

We plot in figure 2 the probabilities P(v® — > isf snrf ) (red curve) and P{v® — > i7© surf ) 
(blue curve) with respect to neutrino energy; where by v® and ^® surf we mean respectively 
the neutrino state at the Sun's center and surface. In this plot we assumed (sf 4 = 0.03, #24 = 
0,fl 34 = 0) for uf -> vf swi ; and (<9 14 = 0,s 2 24 = 0.01,034 = 0) for v® -> P s 0surf ; and 
Am^j = 1 eV 2 for the both curves. As can be seen from the curves, the onset of the active- 
sterile flavor conversion occurs at the energy given by eq. (2.3), that is ~ 85 GeV for v e and 
~ 240 GeV for i T . The rapid oscillation in both of the curves in figure 2 are due to the 1-4 
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^v(GeV) £y(GeV) 



Figure 3. The probabilities P{v® -> ^ 0surf ) (left) and P{v® -> P,f surf ) (right) . For the left plot we 
assumed (sf 4 = 0.03, s| 4 = s| 4 = 0) and for the right plot (s 2 4 = 0, = 0.01, s| 4 = 0); and for the 
both plots Am 41 = 1 cV 2 . 



mixing. The vacuum oscillation length of 1-4 mixing is given by 

~ Am^ ~ 246 ^ ( M6ESV J ( A^j ' <2 ' 4) 

The 14-induced oscillation contribute the two following terms: sin 2 2(9i 4 sin 2 (7ri?0/LQg C ) and 
sin 2 26*24 sin 2 (7ri? Q /LQg C ) to the f® — > i/ 0surf and P® — > z/ Qsurf oscillations, respectively; where 
Rq ~ 7 x 10 5 km is the radius of Sun. In figure 2 we assumed sin 2 26>i4 = 0.11 and sin 2 26> 2 4 = 
0.04; which justifies the relative smaller amplitude of rapid oscillation in the P® — > z?© surf 
curve with respect to the vf — > ^ Qsurf curve. 

For the case of electron neutrino at the center of Sun, i/®, with energy E v > 85 GeV, the 
resonant flavor oscillation due to Am^ plays role such that the state of the neutrino at the 
surface of Sun is almost uf smi . The left plot in fi gure 3 shows the numerical calculation of 
probabilities P{y® — > i/® surf ) with respect to the neutrino energy for the mixing parameters 
the same as for red curve in figure 2. The right plot in figure 3 shows the probabilities for 
the case of v® at the center of Sun and again with the same mixing parameters as the blue 
curve in figure 2. For the v® also the MSW resonance makes the state of neutrino almost 
-©surf £ Qr ener gj es > 240 GeV, as can be seen from the plot. However it should be noticed 
that the complete conversion z7® — > z/® surf occurs when S34 = 0. Assuming S34 7^ 0, the state 
of the neutrino at the Sun's surface is the superposition of z/® surf and z^® surf . 

From the surface of Sun to the surface of Earth neutrinos evolve in vacuum. For a 
continuous energy spectrum which is the case for secondary neutrinos produced in the an- 
nihilation modes XX ~^(qQ^ fermions, gauge bosons), all the oscillatory terms between Sun 
and Earth average out and the oscillation probabilities are the same as at the surface of 
Sun. However, for the monochromatic neutrinos from the annihilation modes XX ~^ v ^ t ne 
oscillatory terms can play role which we will discuss in the next section. 

An important feature of this resonant conversion is that it take place even for very 
small values of mixing angle #14, #24 and #34. Considering the adiabaticity of conversion, it 
can be shown that resonance occur for (s 2 4 , s| 4 , s| 4 ) > 10~ 3 , which is far below the present 
upper limits from the short baseline experiments. Thus, we propose that in the indirect DM 
searches in the neutrino telescopes such as IceCube, this point should be considered. In the 
next section we discuss this point for the case of monochromatic neutrinos in more detail. 
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3 Oscillation probability at the Earth for monochromatic neutrinos 



m x and thus 



The neutrinos produced in XX ~^ y v are monochromatic with the energy E, 
the spectrum of neutrinos at the production point is a sharp line. It has been shown in 
[14] that thermal velocity distribution of x> natural wave-packet width of neutrinos at the 
production point and decoherence due to gravitational acceleration of x do not smear the 
sharp line feature of the spectrum. The other two potential sources of smearing are neutral 
current interaction of neutrinos with the nuclei in Sun's medium and the r-regeneration. 
The neutral current interaction with nuclei N in the Sun's medium could change the energy 
of neutrinos and so smears the sharp line in the spectrum; but since the cross section of 
neutral current scattering u{k)N — > v(k')X is finite in the limit k' — > k (zero momentum 
transfer or forward scattering), the final energy spectrum of neutrinos would be a sharp 
line superimposed on a tail corresponding to scattered neutrinos. The r-regeneration effect, 
which is the reproduction of v e (P e ), [p^) and v r (u T ) from the prompt decay of t + (t~) 
produced in the charged current interaction of v T (y T \ also smears the sharp line in the 
spectrum. However, in the low energy range (m x < 500 GeV) these effects are negligible and 
we focus just on the sharp line in the spectrum. A complete analysis including these effects 
and also annihilation modes XX bb, W + W~ , . . . will appear in [15]. 

In this section we consider the flavor content of the neutrino beam arrived at the Earth 
from the annihilation modes XX ~ > u a^a where (a = e, /x,t) . Also we consider the anni- 
hilation XX ~~ u sV s which is a new possibility in the presence of sterile neutrinos. Basically, 



the amplitude A of flavor oscillation v® 
at the Earth is given by 



,Det 



from the center of Sun to the detector place 



M», 







Dot 



jf surf )exp 



2E„ 



M»\ 



5surf 



(3.1) 



where f® surf represent the neutrino mass eigenstates at the surface of Earth and v^ ct is the 
neutrino flavor state at the site of IceCub detector at South Pole. The exponential factor 
comes from the mass eigenstate's evolution in vacuum between the surface of Sun and the 
surface of Earth; and Les is the Earth-Sun distance ~ 1.5 x 10 8 km. The probability of flavor 



oscillation is given by \A{v® — > )| . Because of the large distance £es> the exponential 
factor due to Am^ average out from the oscillation probability; but, as it shown in [8, 14], 
the exponential factors due to 12-mixing and 13-mixing do not wash out. 

Let us first consider the three active standard neutrino oscillation. It should be noticed 
that unlike the familiar solar neutrinos from nuclear fusion in the Sun, with the energy 
E v ~ MeV, the neutrinos from \X annihilation with energy E v ~ 100 GeV do not leave the 
surface of the Sun in mass eigenstates. In the case of solar neutrinos from nuclear fusion, the 
state of the neutrino at surface of Sun is practically v%\ but, in the case of neutrinos from XX 
annihilation, the state of the neutrino at the surface of Sun is a coherent superposition of v\ 
and V2 (assuming #13 = 0). Thus, generally in the case of monochromatic neutrinos (which 
is the case for XX annihilation directly to neutrinos) the oscillatory terms in eq. (3.1) do not 
average out. It was shown in [8, 14] that this non-averaging of oscillatory terms leads to a 
seasonal variation of the flavor oscillation probability due to the change in the Earth-Sun 
distance during the year (which is ~ 5 x 10 6 km). 



Recently the authors of [21] proposed a new method to distinguish between these annihilation modes based 
on the r-regenerated neutrinos. 
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21 Mar 1 May 1 June 1 July 1 Aug 1 Sep 21 Mar 1 May 1 June 1 July 1 Aug 1 Sep 



Time of year Time of year 

(a) (b) 

Figure 4. Seasonal variation of the probability P{vf —> v^ ct ) (a) and P{y® — > v^ ot ) (b). For figure 
(a) and (b) we assumed (sf 4 = 0.1, s\ A = s§ 4 = 0) and (s 2 4 = s| 4 = 0.1, s§ 4 = 0), respectively; and 
for both figures: Am^ = 1 eV 2 and E v = 300 GcV. 

In 3+1 scenario, the MSW resonance discussed in previous section would change the 
seasonal variation. Let us temporarily ignore the present upper bounds on the active-sterile 
mixing parameters for illustrative reasons. For example, consider the case of \X ~~ v e^e with 
m x > 85 GeV. The v ® state evolving toward the surface of the Sun undergo resonance and 
the state of neutrino leaving the surface of Sun is mostly vf snri . Thus, the sum in eq. (3.1) 
reduces to one term and taking the absolute value to obtain the probability wash out the 
exponential term. However, inside the Earth, neutrinos propagate different paths toward 
the IceCube site during the year due to the change in Sun's zenith angle with time. We 
consider the period of the southern fall and winter seasons (21 March-23 September) where 
the Sun is below the horizon at the IceCube site and thus the background from atmospheric 
muons is suppressed by the Earth's matter. The zenith angle of Sun is 90° at 21st of March 
and takes its maximum value at 21st June (winter solstice in southern hemisphere) which 
is ~ 113.2°. As a result, the propagation length of neutrinos inside the Earth (from the 
surface of Earth at Sun's direction to the IceCube site) is zero at 21st March and 23rd Sep 
and take its maximum value ~ 5100km at 21st June; and considering the L^ c in eq. (2.4), 
we expect to observe the oscillation pattern of 1-4 mixing at the IceCube location. Thus, 
for example, in the case of P{y® — > f,P et ) we expect a constant probability with respect to 
time which is modulated by the 14-induced oscillation inside the Earth. To illustrate this 
argument, we plot in figure 4(a) the probability P(yf — > v^* ct ) with respect to time obtained 
by numerical solution of the evolution equations from the Sun's center to the IceCube site 
at South Pole. For this figure we assumed (sf 4 = 0.1,024 = 0,^34 = 0), Am^ = leV 2 and 
E u = 300 GeV. As can be seen, the pattern of the probability is a constant (with respect to 
time) curve modulated by an oscillatory component. The constant value is sf 4 = 0.12 which 
comes from the decomposition of f® surf to ^® surf and the oscillatory pattern comes from 
the 14-induces oscillation inside the Earth with the amplitude 4cf 4 s^ 4 ~ 0.05. The nearly 
plateau region in the middle of figure 4(a) comes from this fact that at the winter solstice in 
southern hemisphere, the zenith angle of Sun viewed from the South Pole reach it maximum 
value and after this date it begins to decrease. Thus, the small rate of change in the Sun's 
zenith angle which leads to a small rate of change in the path length of neutrino propagation 
inside the Earth, suppresses the oscillation pattern at solstice dates. 
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Figure 5. Seasonal variations of P(v® — > v^ et ) (a) and P{v® — > v^ et ) (b). In each plot the red (solid) 
curve is for 3v standard oscillation and the blue (dashed) curve is for 3+1 scenario. In both plots we 
assumed E v = 300 GeV and for the 3+1 scenario parameters we assumed: s\ A = s| 4 = s 2 ^ = 0.008 
and Ami, = 1 eV 2 . 



Let us now consider the probability P(v® —> vf? et )- In this case again the state of the 
neutrino at the surface of the Sun is almost i/® sur and from eq. (3.1) we do not expect any 
seasonal variation in the probability. However, at the surface of the Earth, the vf SUTi state 
can be decomposed to flavor states z/® surf , i/® surf and f® surf with coefficients sn, C14S24 an d 
C14C24 respectively (assuming #34 = 0). It should be noticed that if S24 = 0, the probability 
P(y® — > ^? et ) would be zero; and because of this we assume S24 7^ in this case. As 
discussed above, due to the comparable oscillation length L^ c with the propagation length of 
neutrinos inside the Earth, oscillatory pattern from j/© surf — > iP ct is expected (the oscillation 
probabilities of z/® surf — y iP ct and i/© surf — y u9 et are very small). Figure 4(b) shows the 
seasonal variation of P{v® — y v^ ct ). In this figure we assumed (s 2 4 = s\ 4 = 0.1, s| 4 = 0), 
Am 2 41 = leV 2 and E v = 300 GeV. In this case the amplitude of oscillation is the product of 
c i4 s 24 ~ 0- 1 (from vf smf decomposition to f® surf ) by the factor 4c| 4 c| 4 s| 4 ~ 0.32 (from the 
oscillation probability induced by 14-mixing inside the Earth), which is smaller in comparison 
with P{vf — > v^ et ) in figure 4(a). 

The plots shown in figures 4(a) and 4(b) are for illustrative purposes and the values 
assumed for the 3+1 mixing parameters are larger than the current upper limits. With 
the assumption of the best-fit values for s 2 4 ,S2 4 and Am 41 from the short baseline neutrino 
oscillation experiments [12] and also the upper limit on s 2 4 from MINOS experiment [19], 
practically the oscillation probabilities P{v® — y v^ et ) and P(p®, — y i>j? ct ) suppress to zero. 
However, it should be noticed that a nonzero value for all the mixing parameters S14, S24, S34 
leads to a nonzero decomposition of z/® surf to all the flavor states, which open new oscillation 
channels u® suri -+ u^ et and j/© surf ->■ z/°et_ 

It should be mentioned that the main detection mode of neutrinos in the energy range 
~ 300 GeV in neutrino telescopes such as IceCube is through /x-track events. The main 
contribution to /u-track events come from the charged current interaction of and with 
the nuclei inside or in the vicinity of IceCube site. Thus, the IceCube is mainly sensitive to 
the and flux at its site. 

We have shown in figures 5(a) and 5(b) the seasonal variation of P{v® — > ^ ct ) and 
P{y® — > z^ et ), respectively. Figures 6 and 7 also show seasonal variation for muon and 
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Figure 6. Seasonal variations of P(v® — > v^ et ) (a) and P{v® — > v^ et ) (b). In each plot the red (solid) 
curve is for 3v standard oscillation and the blue (dashed) curve is for 3+1 scenario. In both plots we 
assumed E v = 300 GeV and for the 3+1 scenario parameters we assumed: s\ A = s| 4 = Sg 4 = 0.008 
and Am?-, = leV 2 . 




21 Mar 1 May 1 June 1 July 1 Aug 1 Sep 21 Mar 1 May 1 June 1 July 1 Aug 1 Sep 

Time of year Time of year 



(a) (b) 



Figure 7. Seasonal variations of P{v® —> v^ et ) (a) and P{v® —> v^ et ) (b). In each plot the red (solid) 
curve is for 3v standard oscillation and the blue (dashed) curve is for 3+1 scenario. In both plots we 
assumed E v — 300 GeV and for the 3+1 scenario parameters we assumed: s 2 4 = = = 0.008 
and Ami, = 1 eV 2 . 



tau (anti)neutrino at the center of Sun, respectively. In each plot the blue (dashed) curve 
represents the seasonal variation in 3+1 scenario and the red (solid) curve corresponds to 
seasonal variation in 2>u standard oscillation; where in both cases E u = 300 GeV. For the 
3+1 scenario in these plots we assumed Am 41 = leV 2 and s\ A = s\± = s^ 4 = 0.008; which is 
well below the current upper limits. As can be seen, for figures 5(a), 6(b) and 7(b), where 
the MSW resonance occur in the Sun, the probability is suppressed; but however as we 
mentioned the nonzero values of s| 4 and s| 4 leads to a small contribution to the probability 
in figures 6(b) and 7(b). The oscillatory patterns in figures 6(a), 7(a), 6(b) and 7(b) are 
the modulation of two components: i) the short wave-length oscillation induced by Awijj 
between the Sun and Earth; ii) the long wave-length oscillation induced by Am|i between 
the Sun and Earth. The oscillation length of 1-2 and 1-3 mixings at E v = 300 GeV are 
~ 10 7 km and ~ 3 x 10 5 km, respectively. Thus, considering the change in Les during the 
southern fall and winter which is ~ 2.5 x 10 6 km, justifies the oscillatory patterns. It should 



- 9 - 



1.0 




o.o L 



2 1 Mar 1 May 1 June 1 July 1 Aug 1 Sep 

Time of year 



Figure 8. Seasonal variation of the probability P(v® — > v^ ct ) (red solid curve) and P{vf — > v^ et ) 
(blue dashed curve) for E v = 400 GeV. For this figure we assumed s\ A — s§4 = s 34 = 0.008 and 
Amlx = lcV 2 . 



be noticed that in the case of figures 6(b) and 7(b), the 1-4 mixing also induces oscillation 
inside the Earth and as can be seen in these figures the blue (dashed) curves, corresponding 
to 3+1 scenario, represent irregularities. 

The other interesting case is the annihilation XX v s v s in the center of Sun (for a model 
predicting this annihilation see [6]). Figure 8 shows the seasonal variation of P(yf — > z^ et ) 
(red solid curve) and P (uf — > ^P et ) (blue dashed curve) for E v = 400 GeV. For this figure 
also we assumed Am^ = 1 eV 2 and s\ A = s\ A = s| 4 = 0.008. For uf created at the center 
of Sun, the state of the neutrino leaving the surface of Sun is almost ci2|^® surf ) + •si2|z- / 2 )sur ) 
and so according to eq. (3.1) the seasonal variation due to Am 2 ^ is expected to be seen. 
The minimum and maximum of the amplitude of this oscillation with respect to time come 
from the decomposition of c 12 \vf snrf ) + s 12 \v^ smi ) and s 12 \v? suii ) + ci 2 |^ surf ) to the state 
|f® surf ), respectively. It is straightforward to see that the maximum (^c^c^cf^s 2 ^ ~ 0.28) 
and minimum (s\ A s\ A ~ 0) matches with the red (solid) curve in figure 8. For the case of 
iff, the neutrino state at the surface of Sun is mostly |^f surf ) and so we expect a constant 
oscillation probability c^^s^ ~ 0.5 (assuming S13 = 0) which is in agreement with the blue 
(dashed) curve of figure 8. The small amplitude oscillation pattern in the blue (dashed) curve 
are due to the modulation of two components: the 14-induced oscillation inside the Earth 
and 13-induced oscillation from the surface of Sun to the surface of Earth which is a result 
of the fact that the state of the neutrino at the surface of the Sun is not exactly \u% suri ) and 
a few percent |^? surf ) admixture exists. 



4 Conclusion 



We have considered the evolution of neutrinos with energy E u ~ 100 GeV inside the Sun in the 
presence of new sterile neutrino state. Neutrinos in this energy range can be produced from 
the annihilation of DM particles (x) gravitationally trapped inside the Sun. Accumulation 
of DM particles in the Sun rises the number density which result in sizable DM annihilation 
rate. Neutrinos are the ubiquitous final product of DM annihilation either with a continuous 
energy spectrum from the annihilation modes XX ~ * QQi W + W~ , ... or with a monochromatic 
spectrum in the direct annihilation XX ~ * 
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The recent global analysis of the data from short baseline neutrino oscillation experi- 
ments [12] favor the presence of one (or more) mostly sterile neutrino mass eigenstates with 
mass splitting Am 2 ~ leV 2 . The new mass splitting Am 2 ~ 1 eV 2 would result in resonant 
flavor conversion inside the Sun. We have shown that the resonant flavor conversion would 
occur for v e with energies E v > 85 GeV and P^/ T with energies E v > 240 GeV. This antic- 
ipated MSW flavor conversion deplete the flux of v e and z^/ T produced in the annihilation 
of DM particles. It should be noticed that the resonance and the subsequent depletion of 
neutrino flux take place even for very small active-sterile mixing angles, far below the present 
upper bounds from short baseline experiments. Thus, we propose that this point should be 
considered in the analysis of data from neutrino telescopes aiming at DM search from Sun. 
Interestingly, even non-observation of the expected excess of events from DM annihilation in 
the Sun's direction in neutrino telescopes can be interpreted by this phenomenon. 

As an example of this effect, we considered the case of monochromatic neutrinos pro- 
duced in the direct annihilation of DM particles to neutrinos XX vv- It was shown in 
[8, 14] that for the monochromatic neutrinos, the oscillation probabilities and as a result the 
number of events in neutrino telescopes exhibit seasonal variation. We calculated the oscilla- 
tion probabilities with the assumption of presence of sterile neutrino state. As we expected, 
the depletion of neutrino flux at the Earth is severe for the v e and z/^/ r at the center of Sun. 
A new annihilation channel which emerge in the presence of sterile neutrinos is XX u sV s - 
We calculated the oscillation probabilities for this annihilation mode at the Earth. In this 
case, again due to the resonant flavor conversion at the Sun, the probabilities v® — > u^ ct 
and vf — > 9^ et would be significant in spite of the strong upper limits on the sterile-active 
mixing angles. 
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